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SUMMARY 


In  many  equipments,  boiling  occurs  in  confined  space  such  as  the 
clearance  between  the  tube  and  the  support  plate  of  steam  generators. 
Corrosive  concentration  builds  up  at  the  boundary  of  dryout  zone  and  induces 
severe  damage.  The  knowledge  on  this  kind  of  boiling  phenomena  is  very 
limited.  It  is  the  objective  of  this  research  to  understand  this  fundamental 
heat  transfer  of  this  problem  through  systematical  analysis  and  experimental 
studies. 

This  report  describes  the  progress  of  research  in  the  past  year  on 
the  following  areas: 

Chapter  2.  Analysis  of  single  phase  forced  convection  in  the  heat  exchanger 
crevices  between  the  tube  and  baffle-plate  using  creeping  flow 
approximation  for  concentric,  eccentric,  and  inclined  annular 
geometries. 

Chapter  3.  Analysis  of  convective  boiling  and  dryout  in  the  support-plate 
crevices  in  steam  generators.  Compared  with  experimental  data, 
and  recommended  for  designs  to  avoid  the  dryout. 

Chapter  4.  Experiment  of  boiling  and  dryout  in  crevices  with  closed  bottom 
for  the  application  of  tube-sheet  crevices.  Boiling  in  confined 
space  is  studied  in  detail  with  various  fluids.  The  dryout  in 
confined  space  is  also  studied  for  various  geometry  and  fluids. 

A  semi -empirical  correlation  Is  established. 

Chapter  5.  The  current  status  of  forced  convective  boiling  and  dryout  in 
crevices  is  presented.  The  obtained  data  are  discussed. 
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CONCLUSION 

The  fluid  flow  and  heat  transfer  of  the  leakage  stream  between  tube  and 
baffle-plate  in  heat  exchangers  are  studied  for  various  possible  geometries  as  ^ 

indicated  in  Chapter  2.  The  information  would  be  helpful  for  better  design  and 
performance  evaluation  of  heat  exchangers  in  general. 

The  boiling  and  dryout  in  support  plate  crevices  can  be  calculated  easily  j 

following  the  computation  scheme  stated  in  Chapter  3.  The  computation  has  been 
checked  with  experimental  data.  Better  design  and  performance  evaluation 
through  computation  becomes  possible.  To  avoid  dryout  completely,  square  hole 
or  holes  with  large  difference  in  curvatures  at  contact  point  are  recommended. 

The  boiling  in  confined  space  with  closed  bottom  generally  has  high 
heat  transfer  coefficient  but  low  critical  heat  flux.  The  Bond  number  is  impor¬ 
tant  in  characterizing  the  boiling  regimes.  At  critical  heat  flux,  the  mechanism 
is  controlled  by  counterflow  flooding.  A  correlation  for  various  geometry  and 
fluids  is  proposed.  In  general,  it  is  independent  of  Bond  number  as  inferred  j 

by  the  flooding  correlation.  Based  upon  this  general  correlation,  the  dryout 
of  steam  generator  crevices  in  the  tube-sheet  can  be  estimated. 

The  convective  boiling  heat  transfer  in  confined  space  has  been  reported  ) 

here.  It  is  similar  to  conventional  boiling;  however,  the  heat  transfer 
coefficient  is  high  and  the  critical  heat  flux  is  low.  Further  study  are 
underway.  J 
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CHAPTER  1.  INTRODUCTION 


Boiling  at  conventional  heat  transfer  surface  has  been  studied 
extensively  in  the  past  thirty  years  [1].  However,  at  the  same  time, 
the  boiling  in  confined  space  has  been  almost  completely  neglected.  The 
existing  information  [2],  [3]  is  limited  and  incomplete. 

In  many  equipments,  boiling  occurs  in  confined  space.  For  example, 
in  the  steam  generators,  clearance  is  allowed  between  the  heated  tube  and 
its  support-plate  where  the  tube  runs  through.  In  the  clearance,  flow  is 
reduced  and  boiling  occurs,  which  is  shown  in  the  Fig.  1(a).  When  the 
liquid  is  highly  subcooled,  single  phase  forced  convection  occurs  in  the 
crevice.  The  behavior  of  the  leakage  stream,  as  shown  in  Fig.  1(b)  is  of 
great  interest  to  the  designer  of  heat  exchangers.  At  the  tube  sheet  of 
steam  generator  a  crevice  with  closed  bottom  is  formed.  The  boiling  and 
dryout  may  also  occur  in  this  tube  sheet  crevice  as  indicated  in  the  Fig.  1(c). 

The  boiling  of  liquid  in  confined  space  leads  to  permanent  dryout. 
Corrosive  concentration  may  build  up  at  the  boundary  of  dryout  zone  where 
boiling  occurs.  Solid  deposits  may  also  occur  at  the  dryout  boundary. 

The  understanding  of  the  thermal -hydraulic  related  corrosion  in  confined 
space  is  lacking.  This  is  primary  because  the  boiling  heat  transfer  in 
confined  space  is  not  really  known.  The  fundamental  understanding  of 
the  boiling  heat  transfer  in  confined  space  is  also  important  to  the  advances 
in  power  engineering  and  lubrication  engineering. 

In  view  of  this  need,  a  systematic  study  of  boiling  heat  transfer  in 
confined  space  is  performed  at  the  Department  of  Mechanical  Engineering 
of  Carnegi e-Mel  Ion  University  under  the  support  of  the  Office  of  Naval 
Research.  Both  experimental  and  analytical  research  are  performed. 
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The  progress  of  the  research  has  been  presented  in  this  report 
in  the  following  chapters: 

Chapter  2.  Analysis  of  Single  Phase  Convection  in  Support-Plate  Crevices. 
Chapter  3.  Analysis  of  Convective  Boiling  and  Dryout  in  Support-Plate 
Crevices. 

Chapter  4.  Experiments  of  Boiling  and  Dryout  in  Annular  Crevices  with 
Closed  Bottoms. 

Chapter  5.  Experiments  of  Convective  Boiling  and  Dryout  in  Annular 
Crevices. 
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CHAPTER  2 

ANALYSIS  OF  SINGLE  PHASE  CONVECTION 
IN  SUPPORT-PLATE  CREVICES 


ABSTRACT 


The  laminar  flow  and  heat  convection  in  annular  type  crevices 
of  various  eccentricities  and  inclinations  are  analyzed.  The  annulus 
is  heated  by  hot  flowing  fluid  at  the  inside  of  the  inner  tube.  The 
flow  field  is  calculated  using  creeping  flow  approximation  of  boundary 
layer  theory.  The  thermal  field  is  studied  two-dimensionally  in 
terms  of  the  local  bulk  mean  temperature  of  the  fluid.  Calculation 
is  performed  numerically  but  in  non-dimensional  form.  It  is  observed 
that  the  angle  of  inclination  influences  the  flow  and  heat  transfer 
significantly  especially  when  the  eccentricity  is  severe.  When  the 
angle  of  inclination  increases,  the  flow  rate  is  reduced,  the  bulk 
mean  temperature  at  the  exit  of  crevice  increases,  and  the  averaged 
heat  transfer  coefficient  in  crevice  decreases.  Comparisons  are 
made  with  respect  to  the  results  of  concentric  annulus. 
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NOMENCLATURE 

A  Cross  section  area  of  axial  flow  in  crevice 
c  Averaged  gap  thickness  in  the  crevice 
Cp  Specific  heat  of  the  fluid  in  crevice 
e  Eccentricity,  defined  in  Figure  1 

f  Apparent  friction  factor  for  crevice,  define  in  equation  (41) 
hj  Heat  transfer  coefficient  inside  the  tube 
hQ  Local  heat  transfer  coefficient  at  outside  of  tube  in  crevice 
F  Averaged  heat  transfer  coefficient  on  the  tube  in  the  whole  crevice 
H  Local  overall  heat  transfer  coefficient  from  the  inside  of  the  tube  to 
crevice  defined  in  equation  (17). 

Thermal  conductivity  of  the  fluid  in  crevice 
ks  Thermal  conductivity  of  the  tube  material 
L  The  axial  length  of  *ube  in  crevice 

m  The  total  mass  flow  rate  through  the  crevice 

n  Index  of  inclination  of  the  tube 

Nu*  Characteristic  Nusselt  number  of  crevice,  defined  in  equation  (38) 

Nu  Averaged  Nusselt  number  in  crevice,  defined  in  equation  (44) 
p  Local  pressure  of  fluid 

pQ  Pressure  difference  across  the  crevice 

P  Non-dimensional  pressure  defined  in  equation  (24) 

Pe*  Characteristic  Peclet  number  of  crevice,  defined  in  equation  (39) 
q^  Averaged  heat  flux  of  the  tube 
R  Radius  of  the  tube 

s  The  thickness  of  the  tube 

S  Ratio  of  the  dimensions  of  the  flattened  crevice,  defined  in  equation  (22) 

T  Local  temperature  of  fluid 
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T  Local  bulk  mean  temperature  of  the  fluid  in  crevice 
T.j  Bulk  mean  temperature  of  hot  fluid  inside  tube 

Tq  Temperature  of  fluid  at  inlet  of  crevice 
Te  Bulk  mean  temperature  of  fluid  at  the  exit  of  crevice 
u  Local  velocity  of  fluid  in  8  direction,  a  function  of  z 

u  Local  velocity  of  fluid  in  8  direction,  averaged  over  the  gap 

U  Non-dimensional  velocity  of  fluid,  in  e  direction  defined  in  equation  (25) 
v  Local  axial  velocity  of  fluid,  a  function  of  z 

v  Local  axial  velocity  of  fluid,  averaged  over  the  gap 

v*  Characteristic  axial  velocity  of  fluid,  defined  in  equation  (40) 

V  Non-dimensional  axial  velocity  of  fluid,  defined  in  equation  (26) 
x  Coordinate  measured  in  e  direction 

X  Non-dimensional  x  coordinate,  defined  in  equation  (20) 
y  Coordinate  measured  in  axial  direction 

Y  Non-dimensional  y  coordinate,  defined  in  equation  (21) 

z  Coordinate  measured  normal  to  tube  surface  across  the' gap 
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Greek 

a  Thermal  diffusivity  of  fluid  in  crevice 
fi  Local  thickness  of  gap  in  crevice 

a  Non-dimensional  gap  thickness  of  crevice,  defined  in  equation  (23) 
e  Eccentricity  ratio  defined  in  equation  (2) 
y  Angle, defined  in  Figure  1 

e  Non-dimensional  local  bulk  mean  temperature  of  fluid,  defined  in 
equation  (27) 

u  Viscosity  of  fluid  in  crevice 
v  Kinematic  viscosity  of  fluid  in  crevice 
p  Density  of  fluid  in  crevice 


INTRODUCTION 


The  fluid  flow  and  heat  transfer  behavior  in  concentric  annulus  are 
well  known  [1].  However,  in  many  practical  thermal  systems  the  axis  of 
the  inner  tube  may  be  inclined  with  respect  to  the  outer  shroud.  Since  the 
annular  gap  could  be  very  narrow,  if  the  inner  tube  is  slightly  off-centered 
the  periferal  variation  of  the  flow  channel  gap-thickness  could  be  very  severe. 
The  fluid  flow  and  heat  transfer  in  this  kind  of  annular- type-crevices  will 
be  relatively  complex.  In  conventional  shel  1-and-tube  heat 

exchangers  and  steam  generators. the  hot  tube  runs  through  the  hole  in  baffle- 
plate  or  support- pi  ate  with  a  narrow  mechanical  clearance  in  between  of  the 
tube  and  wall  of  the  hole.  Due  to  the  misalignment  of  the  plates  and  the 
displacement  of  tube  in  the  flow  field  it  is  likely  that  the  tube  may  be 
inclined  in  the  hole  to  form  an  inclined  annular-type  crevice. 

The  flow  is  forced  through  the  crevice  by  the  pressure  difference 
across  the  plate.  The  fluid  is  also  heated  by  the  hot  tube  in  the  crevice. 

The  knowledge  of  the  flow  rate  and  temperature  of  the  leakage  stream  is 
of  great  importance  to  the  understanding  of  the  thermal  performance  of  these 
'heat  exchange  equipments.  This  type  of  heat  transfer  problems  also  occur  in 
many  other  industrial  applications. 

The  analysis  of  the  fluid  flow  in  concentric  annular  orifices  has  been 
reported  in  references  [2,3J.  The  fluid  mechanic  study  of  annuli  of  small 
clearances  has  also  been  performed  in  references  [4,5]  for  concentric 
and  eccentric  configurations.  However,  no  analysis  has  been  conducted  for 
the  general  case  of  annular-type  crevices  with  the  inner  tube  inclined  and 
placed  eccentrically.  The  analysis  of  the  flow  in  this  generalized  condition 
is  one  of  the  objectives  of  the  present  study. 
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The  other  objective  of  the  present  study  is  to  analyze  generally 
the  convective  heat  transfer  in  annular-type  crevices.  In  most  applications 
the  tube  is  heated  from  inside  by  another  flowing  hot  fluid  such  that  the 
outside  surface  condition  of  the  tube  is  neither  at  a  constant  heat  flux 
nor  at  a  constant  temperature.  The  heat  transfer  of  eccentric  annulus 
at  constant  heat  flux  has  been  studied  in  I6j  for  slug  flow.  However,  no' 
analysis  has  been  performed  for  laminar  flow  in  annular-type  crevices  with 
inclined  inner  tube.  In  fact,  this  information  is  of  great  importance  to 
the  heat  exchanger  designs  [7J. 

In  the  following  sections  the  fluid  flow  and  heat  convection  of  annular- 
type-crevices  will  be  analyzed.  To  limit  the  complexity  of  the  problem  the 
following  assumptions  were  made: 

Cl)  The  problem  is  steady  state  with  constant  properties  for  the  fluid. 

(2)  The  flow  is  laminar  in  the  crevice  with  negligible  inertial 
•  effect  and  negligible  energy  dissipation.  The  entrance  effects 
are  neglected. 

C3)  Heating  is  provided  by  the  hot  fluid  in  the  tube.  The  tube  is 
generally  not  in  contact  with  the  shroud.  Should  contact  occur, 
the  contact  will  be  limited  to  a  point. 

FORMULATION 

Geometry:  To  describe  the  geometry  of  the  annular- type-crevices 
accurately  the  local  gap  thickness  of  the  flow  channel  has  to  be  quantified. 
At  any  cross-sectional  cut  view  of  the  crevice  the  configuration  of  the  flow 
channel  can  be  described  as  an  eccentric  annulus.  For  two  eccentric  circles 
with  eccentricity  e  and  with  their  radii  having  a  small  difference  c,  the  gap 
thickness  5  can  be  evaluated  accurately  but  not  exactly  from  Figure  1  as 
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5tf)  3  C(1-e  COS  y)  (1) 

where  e  =  e/c  (2) 

is  the  eccentricity  ratio  which  varies  between  0  and  1. 

If  the  inner  tube  is  inclined  with  respect  to  the  outer  shroud,  the 
eccentricity  ratio  e  would  be  a  function  of  the  axial  location.  The 
displacement  of  the  inner  tube  can  be  decomposed  as  two  steps.  First,  the 
tube  moves  from  the  concentric  position  to  an  eccentric  position.  Then  the 
tube  is  inclined  with  an  angle.  In  the  present  study,  a  simple  but  typical 
condition  is  considered  that  the  plane  of  inclination  of  the  inner  tube 
passes  through  the  center  line  of  the  outer  shroud.  This  is  shown  in 

i 

Figure  2.  The  general  formulation  for  the  gap  of  crevice  is 

6  =  c{l-e[l-  2n( 1- (y/L) )]  cos  y>  (3) 

where  the  eccentricity  ratio  e  is  the  value  at  the  cross-section  where  } 

y  =  L.  ‘The  n,  called  the  index  of  inclination, controls  the  angle  of 
inclination  of  the  inner  tube.  When  n  equals  zero  the  inner  tube  is  in 
parallel  with  the  outer  shroud.  At  this  moment  the  equation  (3)  is  the  same 
as  equation  (2).  When  the  n  increases  but  with  the  e  fixed,  the  anlge  of 
inclination  of  the  tube  increases  while  its  relative  position  at  y  *  L  is 

A 

not  changed.  When  n  equals  unity,  the  inner  tube  inclines  symmetrically 
in  the  crevice  from  a  side  view.  For  the  extreme  condition  of  e  *  1 
and  n  *  1,  the  inner  tube  contacts  the  outer  shroud  at  two  points  y  -  0 

> 

and  y  »  L  respectively. 

The  cross-sectional  view  of  the  crevice,  as  shown  in  Figure  1,  is 
symmetric  with  respect  to  the  line  A-B.  Therefore,  only  half  of  the 
total  crevice  has  to  be  calculated.  In  the  present  study  the  gap  thickness 
6  is  generally  much  smaller  than  the  tube  radius  R  such  that  the  curvature 
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of  the  crevice  can  be  neglected.  Therefore,  the  crevice  can  be  analyzed 
as  a  flattened  channel  with  the  varying  height  6  as  shown  schematically  in 
Figure  3.  The  dimension  of  the  flat  channel  is  L  in  the  y  direction  and  ttR 
in  the  x  direction  with  a  local  height  S(xy).  The  symmetry  conditions  will 
be  applied  to  the  lines  of  x  =  0,  and  x  *  ttR. 

Flow  Field:  The  creeping  flow  approximation  to  the  boundary  layer  theory  will 
be  applied  to  the  annular-type  crevices.  For  the  laminar  flow  in  narrow  flow 
passage,  the  flow  field  will  be  dominated  by  the  viscous  force  and  the  pressure 
driving-force.  The  fluid  inertia  force  may  be  neglected  [8]  if 


y^L  < 


c  2 

([-)  «1 


(4) 


where  the  v*  is  the  characteristic  velocity  of  fluid  in  the  crevice  and  will 
be  defined  later.  At  this  condition  the  velocity  profiles  become  parabolic. 
The  parabolic  velocity  profiles  can  be  substituted  into  the  equation  of 
continuity  and  integrated  cross  the  gap  from  0  to  6.  The  result  is 


in?  <{3  H>  *  £  ({3  ■  0  (5) 


which  is  the  Reynold's  equation  of  lubrication  [8]  when  the  walls  are  not 
in  motion.  In  this  equation  the  5  is  given  from  the  known  geometry  of 
crevice  and  the  pressure  field  is  unknown.  The  boundary  conditions  are 


P(x,  0)  =  pQ  (6) 

p(x,  L)  =  0  (7) 

where  the  pQ  is  the  driving  pressure  over  the  crevice,  and 

|f  (0,  y)  =  0  (8) 

ff  (itR.  y)  =  0  (9) 


for  the  symmetry  condition  as  described  previously. 
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■  al  Geometry  of  the  Crevice  Considered  in  Analysis 


fig.  3. 


I 


Heat  Transfer:  Considering  the  constant  properties  and  neglecting  the 
viscous  dissipation,  the  energy  equation  of  the  laminar  flow  in  the  crevice 


becomes 


3T  .  „  3T  _  „  ,32T 

-“<^7 


32T  .  32T \ 

7  7 


The  order  of  magnitude  analysis  can  be  performed  to  justify  the  importance 
of  each  term.  Normalize  the  velocity  by  v*,  x  and  y  by  L,  and  z  by  6. 

The  convective  terms  are  of  the  order  1;  the  conduction  in  x  and  y  directions 
are  of  the  order  (v*  L/a)-1;  and  the  conduction  in  z  direction  becomes  order 
of  [(v*L/a) (6/L)  ]  .  Considering  the  criterion  in  equation  (4)  and  assuming 

the  Prandtl  number  of  the  order  1,  the  conduction  in  z  direction  would  be 
dominant.  The  conduction  in  x  or  y  direction  would  be 

negligible  because  the  characteristic  Peclet  number  (v*l/a)  is  generally 
much  larger  than  one.  The  omitting  of  the  conduction  in  x  and  y  direction  is 
further  supported  by  the  particular  type  of  heating  condition  in  the  present 
study.  The  temperature  of  the  hottest  point  in  the  crevice  is  limited  by  the 
hot  fluid  temperature  in  the  tube;  therefore,  a  severe  hot  spot  is  not  likely 
to  occur  and  the  conduction  in  x  or  y  direction  will  be  limited  too. 

The  remaining  energy  equation  can  be  integrated  with  respect  to  z  from  0 
to  <S  and  gives 


J 


1 


) 
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(13) 


Following  the  derivation  in  the  Appendix,  the  first  term  of  the  equation  (13) 
can  be  written  as 


J  “S  d2  ...I-f  If  <T-T!dz 

0  -^o 


(14) 


where  the  T  is  the  local  bulk  mean  temperature  of  the  fluid.  Since  the  gap 
of  the  crevice  is  narrow,  the  local  temperature  of  fluid  is  very  close  to 
the  local  bulk  mean  temperature.  When  the  gap  thickness  6  does  not  vary 
strongly  in  the  crevice,  the  value  of  3u/3x  will  also  be  small.  Therefore 
the  last  term  in  the  equation  (14)  can  be  neglected  in  general.  Approximate 
the  second  term  in  equation  (13)  in  the  same  manner,  the  energy  equation  becomes 

“{lr+  7e§’°£]  <15> 

0 


If  the  crevice  is  heated  from  the  inner  tube  and  insulated  at  the  outer 
tube,  the  last  term  in  the  equation  (15)  can  be  written  as  qw/pCp.  The  wall 
heat  flux  qw  can  be  described  as  H(Ti  -  T).  Therefore  the  energy  equation  is 

-v  <>6> 

where  the  H  is  the  overall  heat  transfer  coefficient  defined  as 


T7  +  fs+T0  <17> 

The  energy  equation  (16)  is  a  first  order  linear  partial  differential  equation 
of  the  hyperbolic  form.  The  initial  conditions  are 
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T 


y=0 


at  inlet,  and 


31 

3x 


x=0 


=  0 


for  the  symmetry. 


(18) 

09) 


Non-Dimensional izatiori:  The  parameters  can  be  non-dimensional ized  using 


and 


X 

”  TT  R 

(20) 

=  2. 

L 

(21) 

II 

U1- 

(22) 

m  6 
c 

(23) 

II 

g\° 

(24) 

_  u  unR 

Po 

(25) 

„  uvL 

2 

(26) 

Po  c 

=  I-Tft 

V  To 

(27) 

Then  the  Reynold’s  equation  and  the  boundary  conditions  become 


/ , 3  3P\  , 
(A  ^)  + 

_L  (A3  1£)  =  o 

3Y  3Y  U 

(28) 

P(X,0) 

=  1 

(29) 

P(X,1) 

=  0 

(30) 

=  0 

(31) 

ffd.Y) 

*  0 

(32) 

with 


The  non-dimensional  velocity  field  is  described  by 


n  -  1  a2  3P 

u "  ■  tia  n 


V  =  - 


1  a2  3 p 

T2A  ST 


The  energy  equation  and  its  initial  conditions  are 

<.2  36  u  30  .  1  ,LwNu*w1  M 
s  U3X  +  v  a?'  I  k){¥i*){]-e) 


with 


0  =  0  at  Y  =  0 
||  =  0  at  X  =  0 


(33) 

(34) 

(35) 

(36) 

(37) 


where 


Nu*  = 

kl 

Pe*  *  *  and 

P  c2 

(38) 

(39) 

o\t 

II 

*> 

(40) 

the  v*  is  the  characteristic  axial  velocity  in  the  crevice  at  laminar  flow. 


NUMERICAL  METHOD 

The  pressure  field  is  calculated  from  the  Reynold's  equation  of  lubrication 
as  shown  in  equations  (28)  to  (32).  To  solve  this  Poisson  equation,  the 
method  of  Alternative  Directional  Implicit  numerical  scheme  [9]  is  used.  The 
implicit  calculations  are  performed  along  horizontal  lines  for  the  whole 
region  of  the  flow  field;  then  the  calculations  are  performed  along  all  the 
vertical  lines  alternatively.  In  the  implicit  scheme  of  the  solution  along 
a  line  the  method  of  Gauss  Elimination  is  used  to  solve  the  system  of  equations. 
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The  thermal  field  is  then  calculated  using  the  known  velocity  field 
which  is  derived  from  the  pressure  field.  In  solving  the  energy  equation 
the  upwind-difference  scheme[10]  has  been  used  to  maintain  a  stable  numerical 
computation. 

The  whole  field  is  typically  discretized  into  20  x  18  meshes.  Approximately 
20  iterations  are  required  in  the  Alternative  Directional  Implicit  calculations 
to  reach  a  converged  solution.  The  result  is  validated  with  hand  calculations. 

The  overall  computer  time  for  the  calculation  is  about  12  cpu  seconds  using  a 
DEC-20  computer. 

RESULTS  AND  DISCUSSION 

To  illustrate  the  result  of  fluid  flow  and  heat  transfer  calculations,  a 
specific  example  is  studied  in  detail.  The  inner  tube  contacts  with  the  outer 
shroud  at  the  top  but  with  a  small  angle  of  inclination.  Fluid  flows  upward 
in  the  crevice.  The  calculated  pressure  field  is  shown  in  Figure  4  in  non- 
dimensional  form.  The  geometric  configuration  of  the  crevice  is  also  specified 
in  this  figure.  The  pressure  field  is  distorted  and  the  flov'  ‘earn  dv.jvts 
from  the  point  of  contact  where  the  gap  thickness  is  the  minimum. 

For  the  same  crevice,  the  thermal  field  is  shown  in  Figure  5  in  non- 
dimensional  form.  Since  the  gap  is  wide  at  the  location  where  the  X  is  large, 
the  flow  is  faster  over  there  and  the  local  bulk  mean  temperature  rises  rela¬ 
tively  slow  along  the  stream.  In  the  region  where  X  is  small,  the  flow  is 
retarded  due  to  the  narrow  gap  thickness.  The  fluid  temperature  rises  fast 
and  reaches  the  temperature  of  hot  fluid  inside  the  tube  in  a  short  distance. 

This  high  temperature  region  is  relatively  large  in  this  example.  If  the 
inner  tube  is  not  in  contact  with  the  outer  shroud,  the  hot  region  would  be 
small  or  inexist. 


o  jr~ 


Contact  Point 
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Many  other  cases  are  studied  for  various  eccentricity  ratio  s  and 
inclination  index  n  for  this  crevice.  The  parameters  of  this  crevice  has 
been  listed  in  Table  1.  For  the  same  crevice,  the  ratio  of  the  mass  flow 
rates  at  various  conditions  and  the  mass  flow  rate  at  concentric  configuration 
are  shown  in  Figure  6  with  the  eccentricity  ratio  e  as  horizontal  coordinate 
and  the  inclination  index  n  as  parameter.  When  n  equals  zero  the  tube  is 
in  parallel  to  the  shroud.  As  indicated  in  Figure  6,  at  a  fixed  e  the  more 
the  inclination  (i.e.,  higher  value  of  n)  the  lower  the  mass  flow  rate.  This 
is  because  the  flow  streams  are  distorted  when  the  tube  is  inclined  in  the 
crevice.  The  extreme  conditions  occur  when  e  equals  unity,  where  the  inner 
tube  contacts  the  outer  shroud  at  the  top  edge.  When  n  equals  zero,  the  tubes 
have  a  line  of  contact  with  shroud  in  crevice.  Large  amounts  of  fluid  passes 
through  the  wide  opening  at  the  other  side  of  the  crevice  and  gives  the 
highest  possible  mass  flow  rate  in  all  the  configurations  (2.5  tiroes  that  of 
concentric  annulus).  If  n  equals  unity,  the  tubes  has  two  contact  points  at 
the  top  and  bottom  of  crevice  respectively.  The  mass  flow  rate  at  this  condi¬ 
tion  is  the  lowest  of  all  the  cases  studied  (about  70  percent  of  the  concentric 
annulus). 

The  overall  flow  resistance  of  the  crevice  can  be  characterized  by  an 
"apparent"  friction  factor  which  is  defined  in  a  conventional  manner  as 

ip*  '4)l  A  <«d 

c  pA* 

It  is  interesting  to  point  out  that  for  a  crevice  with  fixed  dimensions  the 
total  cross-sectional  area  A  of  the  flow  is  a  constant  value  irrespective 
to  the  eccentricity- of  the  tube  in  the  crevice.  The  apparent  friction  factor  of 
the  crevice  is  shown  in  Figure  7  for  various  e  and  n.  The  Reynolds  number  here 
is  based  upon  the  area  averaged  axial  velocity  of  the  fluid  in  crevice.  The  nature 


of  this  figure  is  similar  to  that  of  Figure  6  due  to  their  interrelationship 

in  equation  (41 ) . 

The  non-dimensional  bulk  mean  temperature  of  the  leakage  stream  leaving 
the  crevice  is  shown  in  Figure  8  for  various  values  of  e  and  n.  The  condition 
of  heating  is  described  in  Table  1.  When  the  tubes  are  concentric  (e=0)  the 
leakage  is  as  hot  as  the  temperature  of  fluid  in  the  tube.  If  the  leakage 
stream  is  more  (for  example, £>0. 5  but  n<0.5)  the  temperature 
of  the  leakage  stream  becomes  lower.  However,  when  the  inner  tube  is  further 
inclined  with  the  n  value  equals  or  larger  than  0.75  the  leakage  stream 
could  be  equal  or  less  than  that  of  concentric  annulus  (see  Figure  6)  the 
temperature  of  leakage  stream  will  be  as  hot  as  the  fluid  temperature  inside  the  tube. 

Based  upon  the  above  information  on  flow  rate  as  well  as  temperature 
of  the  leakage  stream,  the  averaged  heat  transfer  from  the  hot  tube  in  the 
crevice  can  be  evaluated,  -hi  may  define  the  averaged  heat  transfer  coefficient 


\  ■  h<Tt  -  V 

By  energy  conservation  for  the  whole  crevice,  the  averaged  heat  flux  is 

_  «C(Te-T,) 


Therefore  the  averaged  Nusselt  number  in  the  crevice  can  be  calculated. 

C„  c  m  W 
Nu  =  nR  L  k. 


The  averaged  Nusselt  numbers  of  crevice  with  various  geometric  configurations 
are  shown  In  Figure  9  as  a  ratio  to  the  Nusselt  number  of  concentric  annulus. 
It  is  interesting  to  notice  that  for  the  condition  of  present  study,  the 
configurations  with  n<0.5  gives  higher  heat  transfer  coefficient  although  the 


Averaged  Nusselt  Number  of  Various  Crevices 
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bulk  mean  temperature  of  the  exit  stream,  as  shown  in  Figure  8,  is  not  as  hot 
as  the  temperature  of  fluid  inside  the  tube.  This  is  mainly  due  to  the  high 
flow  rate  at  these  conditions  as  indicated  in  Figure  6.  The  Nu  is  essentially 
related  with  the  product  of  mass  flow  rate  and  exit  bulk  mean  temperature 
as  stated  in  the  equation  (44). 

CONCLUSION 

The  laminar  flow  and  heat  transfer  in  annular-type  crevices  of  various 
geometric  configurations  are  studied  in  detail.  The  flow  field  can  be 
calculated  using  creeping  flow  approximation.  The  thermal  field  can  be 
formulated  in  two  dimension  based  upon  the  local  bulk  mean  temperature  in 
crevice.  Comparing  with  the  behavior  of  concentric  annulus,  the  eccentric 
annulus  with  little  inclination  of  the  inner  tube  gives  higher  mass  flow  rate, 
lower  bulk  mean  temperature  at  exit,  but  higher  averaged  heat  transfer  coeffi¬ 
cient  in  crevice.  If  the  angle  of  inclination  is  large,  the  behavior  is 
reversed.’  It  is  suggested  that  the  fluid  flow  and  heat  transfer  behavior  of 
confined  geometry  with  other  types  of  configurations  can  also  be  analyzed 
in  a  similar  manner. 
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CHAPTER  3 

ANALYSIS  OF  CONVECTIVE  BOILING  AND 
DRYOUT  IN  SUPPORT-PLATE  CREVICES 


abstract 

Two-dimensional  analysis  is  performed  for  the  two-phase  thermal-hydraulics  in 
annular-type  clearances  between  the  heated  tube  and  the  support  plate  in  steam 
generators  assuming  the  fluid  is  homogeneous,  non-slip,  and  with  uniform  axial 
pressure  gradient.  Numerical  computation  has  been  performed  for  this  pressure¬ 
boundary  problem,  using  an  under-relaxation  method  for  initial  iterations  to  suppress 
passible  numerical  instability. 

Annular-type  clearances  with  various  eccentricities  and  inclinations  are  studied.  The 
boiling  and  dryout  pattens  at  various  conditions  are  presented  and  compared  with 
available  experimental  data.  The  effect  of  the  difference  in  radii  of  curvatures  near 
the  contact  location  to  the  dryout  is  studied.  Recommendation  for  better  design  to 
avoid  excessive  dryout  and  corrosion  is  proposed. 


NOMENCLATURE 


the  average  gap  thickness  of  the  crevice 
hydraulic  diameter,  equals  2h  for  parallel  plates 
friction  factor  of  the  liquid  at  the  same  mass  flow  rate 
local  cross  mass  flux 

local  axial  mass  flow  rate  per  unit  flow  cross  section  area,  this  is 
local  axial  mass  flux 

r 

local  gap  thickness  of  crevice 

heat  transfer  coefficient  in  the  crevice 

heat  transfer  coefficient  in  the  tube 

local  enthalpy  of  fluid 

latent  heat  of  evaporation 

enthalpy  of  saturated  liquid 

total  mass  flow  rate  in  y  direction  in  the  whole  crevice 

inclination  index 

pressure 

total  pressure  difference 
wall  heat  flux  into  the  crevice 

ratio  of  cross  flow  to  axial  flow,  based  on  mass  flow 
radius  of  the  tube 

local  temperature  of  .fluid,  average  over  gap  thickness 

bulk  mean  temperature  in  the  tube 

coordinate 

coordinate  in  axial  direction. 


Greek  Symbols 


local  fluicj  aualitv 


I 


3  thickness  of  the  tube 

(  eccentricity,  defined  after  equation  (1) 

0  angle,  as  shown  in  Figure  2 

p  density  of  the  fluid  (single  phase  or  two  phase) 

/>(  density  of  liquid 

♦  two  phase  friction  multiplier  based  on  f 

INTRODUCTION 

In  many  tube-and-shell  steam  generators  the  heated  tubes  are  fixed  in  space  by 
support-plates.  As  shown  in  Figure  1,  holes  are  drilled  in  the  support-plate  to  allow 
the  tubes  to  pass  through.  Mechanical  clearance  is  provided  between  the  tube  and 
the  support-plate.  Due  to  the  narrow  passage  of  the  clearance,  fluid  flows  slowly  in 
it.  When  the  tube  temperature  is  higher  than  the  saturation  temperature  of  the  fluid 
at  the  shell-side,  boiling  and  possibly  permanent  dryout  may  occur  in  the  crevices. 
Corrosive  concentrations  build  up  at  the  boundary  of  the  dryout  zone  where  solid 
deposits  may  also  appear.  In  fact,  corrosion  and  denting  of  tubes  has  been 
observed  frequently  at  the  support-plate  crevices  in  many  steam  generators. 

Analysis  has  been  performed  to  predict  the  dryout  pattern  in  crevices  [1]; 
however,  the  analysis  is  limited  to  one-dimensional  flow  where  no  cross-flow  is 
allowed  in  the  azimuthal  direction.  Practically,  when  the  axis  of  the  tube  is  not  in 
parallel  to  the  axis  of  the  hole  on  the  support  plate,  the  flow  in  the  crevice  will 
have  an  azimuthal  component.  In  some  cases,  the  tube  may  contact  the  wall  of  the 
support-plate  hole.  Near  the  contact  location  flow  is  greatly  retarded,  and  boiling 
occurs  easily.  The  vapor  generated  from  the  local  boiling  will  also  have  a  tendency 
to  push  the  adjacent  fluid  away  in'  the  azimuthal  direction  and  induces  cross-flow. 
For  a  better  modeling,  two-dimensional  calculation  is  desired.  In  this  paper  the 
boiling  and  dryout  in  crevices  are  calculated  following  an  approach  of  subchannel 
analysis  which  has  been  used  widely  in  tube  bundle  thermal  hydraulic  analysis. 

The  knowledge  of  thermal-hydraulic  phenomena  in  tube  bundles  has  been  greatly 
advanced  in  the  last  decade.  A  number  of  subchannel  computer  programs  have  been 
developed,  for  example,  in  references  [2-3].  Invariably,  most  of  these  methods 
require  the  evaluation  of  empirical  constants,  usually  by  comparison  with 
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experimental  data,  before  quantitative  predictions  can  be  made.  The  considered 
transport  mechanisms  consist  of  turbulent  mixing,  thermal  conduction,  diversion  cross 
flow,  etc.  Attempts  have  been  made  to  use  subchannel  programs  to  evaluate  the 
dryout  in  crevices.  Numerical  instability  prohibits  successful  calculation. 

In  view  of  the  particular  feature  of  the  crevices,  some  simplifications  can  be  made 
to  the  subchannel  calculation.  The  narrow  crevice  suppresses  the  communication 

between  subchannels.  The  turbulent  mixing  and  the  thermal  conduction  among 
subchannels  can  be  neglected  with  respect  to  the  diversion  cross  flow  and  the  heat 
transfer  at  wall  respectively.  The  use  of  empirical  constants  in  calculation  is  also 
undesirable  due  to  lack  of  knowledge  on  the  physical  behavior  in  crevices. 
Therefore,  an  appr  ximate  but  effective  calculation  method  is  proposed  based  upon 
the  principles  whi~h  have  been  used  in  the  SIMPLE  II  Program  [4].  Simplifications 
have  been  made  to  this  computer  program  but  with  the  extension  to  two  phase  heat 
transfer.  Other  assumptions  and  the  derivations  of  the  analysis  will  be  presented  in 
a  latter  section.  The  method  used  to  suppress  the  numerical  instability  will  also  be 
discussed  later. 

GEOMETRY 

0 

To  describe  the  geometry  of  the  annular-type  crevices  accurately  the  local  gap 

thickness  of  the  flow  channel  has  to  be  quantified.  At  any  cross-sectional  cut  view 
of  the  crevice  the  configuration  of  the  flow  channel  can  be  described  as  an  eccentric 
annulus.  For  two  eccentric  circles  with  eccentricity  e  and  with  their  radii  having  a 
small  difference  c,  the  thickness  5  can  be  evaluated  accurately  but  not  exactly  from 
Figure  2  as 

h<0)  *  c(l-«  cos#)  ,  where  (1) 

e  •  el C  (2) 

is  the  eccentricity  ratio  which  varies  between  0  and  1. 

If  the  inner  tube  is  inclined  with  respect  to  the  outer  shroud,  the  eccentricity  ratio 
•  would  be  a  function  of  the  axial  location.  The  displacement  of  the  inner  tube  can 
be  decomposed  as  two  steps.  First,  the  tube  moves  from  the  concentric  position  to 

an  eccentric  position.  Then  the  tube  is  inclined  with  an  angle.  In  the  present  study, 

a  simple  but  typical  condition  is  considered  that  the  plane  of  inclination  of  the  inner 
tube  passes  through  the  center  line  of  the  outer  shroud.  This  i*  shown  in  Figure  3. 
The  general  formulation  for  the  gap  of  crevice  is 


-  42  - 


h  «  c{l-r[l  -  2nd cos  9)  (3) 

where  the  eccentricity  ratio  <  is  the  value  at  the  cross-section  where  y=  L.  The  n, 
called  the  index  of  inclination,  controls  the  angle  of  inclination  of  the  inner  tube. 

When  n  equals  zero  the  inner  tube  is  in  parallel  with  the  outer  shroud.  At  this 
moment  the  equation  (3)  is  the  same  as  equation  (2).  When  the  n  increases  but  with 
the  *  fixed,  the  angle  of  inclination  of  the  tube  increases  while  its  relative  position 
at  y  =  L  is  not  changed.  When  n  equals  unity,  the  inner  tube  inclines  symmetrically 
in  the  crevice  from  a  side  view.  For  the  extreme  condition  of  *  -  1  ano  "  =  i,  the 
inner  tube  contacts  the  outer  shroud  at  two  points  y  *  0  and  y  *  L  respect 

The  cross-sectional  view  of  the  crevice,  as  shown  in  Figure  2,  is  symmetric  with 
respect  to  the  line  A-B.  Therefore,  only  half  of  the  total  crevice  has  to  be  calculated. 

In  the  present  study  the  gap  thickness  h  is  generally  much  smaller  than  the  tube 
radius  R  such  that  the  curvature  of  the  crevice  can  be  neglected.  Therefore,  the 
crevice  can  be  analyzed  as  a  flattened  channel  with  the  varying  height  h  as  shown 
schematically  in  Figure  4.  The  dimension  of  the  flat  channel  is  L  in  the  y  direction 
and  it R  in  the  x  direction  with  a  local  height  h(xy).  The  symmetry  conditions  will  be 
applied  to  the  lines  of  x  «  0,  and  x  *  irR. 

FORMULATION 

Viewing  the  flow  field  from  the  x-,  lane  as  many  individual  flow  subchannels 
along  the  y-axis,  the  conservation  equations  can  be  formulated  for  each  node  of  each  i 

subchannel.  A  schematic  of  the  divided  subchannels  and  nodes  is  shown  in  Figure  5. 
compared  with  conventional  differential  analysis,  the  present  approach  allows  for 
more  simplifications  for  an  efficient  calculation  yet  with  reasonable  accuracy.  Due 
to  the  simplicity,  the  subchannel  analysis  is  able  to  calculate  the  complicated  two  J 

phase  convective  boiling  with  minimal  extra  effort. 

In  the  analysis,  calculation  is  performed  two-dimensionally  on  the  x-y  plane  with  a 
variable  gap  thickness  h.  The  velocity  and  the  temperature  appearing  in  the  analysis  } 

are  the  averaged  value  and  the  bulk  mean  value  respectively  over  the  thickness  h.  It 
is  assumed  that  the  axial  pressure  gradient  is  identical  on  all  the  nodes  which  are  at 
a  same  axial  elevation  [4],  This  assumption  is  reasonable  because  the  magnitude  of 
the  cross  flow  is  usually  smaller  than  the  axial  flow.  In  the  two  phase  region,  the 
fluid  is  assumed  to  be  homogeneous  with  no  slip  between  the  phases. 


Two  phase  flow  is  analyzed  in  the  confined  space.  The  surface  tension  effect  at 
the  interface  is  neglected.  It  is  true  that  the  gap  thickness  is  extremely  small  near 
the  location  of  contact  such  tnat  the  surface  tension  effect  a! h  would  be  significant 
there.  However,  it  is  also  likely  that  dryout  may  occur  mostly  at  the  contact  area. 
Therefore  the  dryout  interface  would  be  at  the  location  away  from  the  surface 
tension  dominated  area  except  at  the  upstream  front  of  the  dryout  zone. 

At  steady  two  phase  flow  the  pressure  force  is  balanced  by  the  frictional  force 
and  the  acceleration  force.  The  gravitational  effect  is  neglected.  As  will  be 
discussed  later,  generally  the  flow  stream  passes  along  the  dryout  boundary  such  that 
the  change  of  flow  quality  is  not  very  severe.  Therefore  the  acceleration  pressure 
drop  is  negligible  compared  with  the  frictional  pressure  drop  which  is  rather 
significant  in  the  narrow  space.  The  pressure  force  is  balanced  by  the  frictional 
force.  That  means 


o- 


Jl  1 

oh  7 


P1 


(4) 


where  the  f^  is  the  friction  factor  for  liquid  flow  with  the  same  mass  flow  rate  as 
the  two  phase  mixture  the  value  of  96/Re  is  used  for  laminar  flow,  D  is  the  local 

h 

hydraulic  diameter  which  equals  2h,  Gy  is  the  local  axial  mass  flux,  and  <!>  is  the  two 
phase  friction  multiplier.  For  homogeneous  non-slip  two  phase  flow,  the  friction 
multiplier  may  be  described  as  [5] 

(5) 
16) 


♦  «  1.0  when  p  <  0 

*  (_1)  when  p  >  0 

P 


where  p  is  the  local  quality  of  the  fluid,  p  is  the  local  density  of  the 
homogeneous  two  phase  fluid. 


At  a  specific  elevation,  the  local  axial  mass  flux  G  in  the  subchannel  (i.j)  can  be 


yi.j 


evaluated  from  equation  (4).  The  total  axial  mass  flow  rste  in  the  crevice  is  the 


summation  of  all  the  axial  flow  rates  over  all  the  i  subchannels 
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m  =  I  G  h  Ax 
^  yi.j  i.j 

Substituting  equation  (4)  into  (7),  the  total  axial  mass  flow  becomes 


■  •  "  lP^  ]  U-  j]  h(  j 
1=1  Ay  1,J  1,J 


Once  the  axial  flows  are  determined,  the  local  cross  flows  can  be  evaluated  from 
the  mass  conservation  at  a  node.  As  shown  in  Figure  5,  the  mass  conservation 
gives 

Gx  h  Ay  =  Gx  Ay  (9) 

i.j  i+i,i  '-1. j 

♦  Gy  h  Ax  -  Gy  h  Ax 
i.j-t  i.j-1  i.j  i.j+i 

where  the  G  and  Gy  are  known  from  the  axial  flow'  calculations,  the  boundary 
.  s  I./ 

condition  of  symmetry  at  x  =  0  gives 

Gx  *  0  at  i  *  1  (10) 

•-  l.i 

All  the  other  cross  flows  can  therefore  be  evaluated. 


Thermal  Field 


The  enthalpy  balance  in  the  control  volume  of  a  node  can  be  made  according  to 
Figure  5  as 


Gy  h  Ax  H  ♦  Gx  h  AyH 

i.j- 1  i.j- 1  i.j- 1  i"  i.j  i~  i.j  i~  i.j 

♦  q  AxAy  *  Gx  h  AyH 

w  i.j  i+ i.j  i.j 

♦  Gy  h  Ax  H 

7  i.J  i.J*1  '.J 


where  G  is  the  local  mass  flux,  H  is  the  local  enthalpy,  and  q  is  the  local  heat 

W 

flux  from  the  wall.  The  signs  used  in  equation  (11)  are  based  upon  the  flow 
directions  indicated  in  Figure  5.  When  cross  flow  is  different  from  those  assumed  in 
Figure  5  the  form  of  energy  balance  equation  will  be  modified  accordingly. 


t 


For  a  steam  generator  the  wall  heat  flux  can  be  evaluated  as 


112) 


qw  *  (T«  •  T,/(Th  *  A*  ~k  ’ 

where  h  is  the  heat  transfer  coefficient,  o  is  the  tube  wall  thickness,  and  the 

subscripts  t  and  s  are  for  the  tube  side  and  the  shell  side  in  crevice  respectively. 

Once  the  local  enthalpy  is  known  the  local  flow  quality  can  be  evaluted  from 

H  *  H  +  B  H  for  0  <  ,?  <  1  '(13) 

l  sat  fj  r 

NUMERICAL  METHOD 

A  computer  program  has  been  established  to  perform  the  above  calculations.  The 
fluid  flow  in  a  crevice  is  driven  by  a  given  pressure  difference  cross  the  support- 
plate.  Iteration  of  the  mass  flow  in  the  crevice  is  necessary  to  match  this  given 
pressure  drop. 

Initially,  the  local  quality  fl  and  the  total  mass  flow  m  are  assumed.  Then 

calculation  is  performed  on  the  nodes  at  a  same  y  elevation.  The  local  friction 

multiplier  ♦  ^  is  evaluated  from  the  assumed  local  quality  and  fluid  properties  using 
equation  (6).  Using  an  assumed  local  friction  factor  f  ,  the  local  pressure  gradient  (- 
Ap/Ay)  at  this  y  elevation  is  calculated  from  equation  (8).  Iteration  of  equation  (8) 
is  necessary  to  obtain  the  converged  f  .  ^'s.  Then  the  nodes  at  the  next  y  elevation 
are  calculated  in  a  similar  manner  to  evaluate  the  local  pressure  gradient  at  this  new 
elevation.  After  calculation  of  the  whole  flow  field,  individual  pressure  drops  at 

each  elevation  are  summed  to  check  against  the  given  pressure  difference  across  the 

support-plate.  Through  comparison,  the  initially  assumed  total  mass  flow  rate  m  is 

updated.  This  iteration  continues  until  the  total  mass  flow  rate  matches  the  given 
pressure  difference  across  the  support-plate. 

The  next  step  is  to  evaluate  the  local  cross  flow  from  equation  (9).  Then  the  local 
enthalpy  is  calculated  from  equation  (11);  and  the  local  flow  quality  is  obtained  using 
equation  (13).  With  the  up-dated  local  fluid  quality  x  at  every  node  the  overall 
calculation  of  flow  and  energy  is  repeated  until  all  the  local  qualities  converge. 

In  the  computation,  numerical  instability  has  been  observed.  When  the  initial  guess 
of  local  quality  is  very  much  away  from  reality,  the  change  of  local  cross  flow 
between  iterations  could  be  significant.  The  variation  of  local  qualities  in  iterations 

also  becomes  large.  Due  to  the  highly  sensitive  response  of  ♦  to  changes  in  quality 
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or  fluid  density,  numerical  instability  occurs.  The  simpliest  way  of  suporessina  the 
instability  is  to  uncer-.-elax  [6]  the  coupling  term  during  the  iteration  between  tne 
energy  equation  and  the  motion  equation.  This  technique  is  found  to  be  effective. 

In  the  numerical  calculations,  the  typical  mesh  is  about  18  nodes  in  x  domain  to 
cover  180  degrees,  and  20  nodes  in  y  direction  to  cover  the  length  of  the  crevice. 
To  allow  the  flow  quality  to  converge  between  the  motion  and  the  energy  equations 
about  30  iterations  are  required.  The  complete  calculation  of  boiling  and  dryout 
patterns  uses  about  24  seconds  of  cpu  time  on  a  DEC-20  computer. 

r 

RESULTS  AND  DISCUSSION 

The  boiling  and  dryout  in  the  crevice  at  line-contact  is  calculated  [7]  for  data 

t 

conditions  similar  to  those  reported  in  reference  [1].  In  that  experiment,  pressurized 
water  containing  impurities  flows  through  the  crevice  for  a  long  duration  of  time. 
The  liquid  in  the  tube  is  at  a  higher  temperature  under  forced  convection.  Boiling 
occurs  at  the  location  where  the  tube  contacts  the  wall  of  the  support-plate  hole. 
Due  to  the  excessive  boiling  and  dryout,  solid  impurities  deposit  at  the  boundary  of 
the  dryout  zone.  The  experimental  condition  which  has  been  used  as  the  input 
information  of  our  calculation  is  listed  in  Table  1. 

Table  1:  INPUT  INFORMATION  FOR  THE  COMPARISON 
OF  CALCULATION  AND  DATA  [1] 


Tube 

I.D. 

1.905 

cm 

Thickness 

1.092 

mm 

Material 

Inconel-600 

Support -Plate 

Hole  Diameter 

1.9634 

cm 

Height 

1.905 

cm 

Shell-Side 

Water.Pressure 

5.442 

MPa 

Inlet  Quality 

0.15 

Pressure  Drop 

2.04 

KPa 

Tube-Side 

Water 

588 

K 

The  calculated  result  and  the  experimental  data  are  compared  in  Figure  6.  The 
result  of  one-dimensional  analysis  performed  in  reference  [1]  is  also  shown  in  the 
figure.  The  result  of  the  present  calculation  is  similar  to  the  result  of  the  one¬ 
dimensional  calculation,  and  both  results  are  reasonably  close  to  the  experimental 
data.  This  is  mainly  because  the  crevice  with  a  line-contact  configuration  does  not 
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DRYOUT  SUBCHANNEL  ANALYSIS 
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Subchannel  Analysis 
One-Dimensional  Analysis  [i] 
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p  =0.15 


x  (degree)  ' 

TV^  6  The  paiierns  or  dryout  and  iiow  quality 
of  line-contact  configuration. 
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have  a  sizable  cross  flow.  (The  tube  is  in  parallel  with  the  hole  axis).  Therefore  the 
difference  between  the  present  analysis  and  the  one-dimensional  analvsis  is  not 

obvious. 

The  behavior  of  dryout  is  also  studied  for  the  condition  of  one-point-  contact  but 
at  various  inclinations.  As  shown  in  Figure  7  the  contact  point  is  at  the  top  of  the 
clearance  at  zero  degree.  With  more  and  more  inclination  (as  denoted  by  higher 
value  of  n)  the  dryout  zone  becomes  less  and  less.  At  the  extreme  condition  of.n  = 
1.0  when  the  tube  contacts  the  shroud  at  two  points  (at  the  top  of  zero  degree  and 
the  bottom  of  180  degree)  another  dryout  zone  also  appears  at  the  location  near  18C 
degrees  with  significant  effects  to  its  downsfeam.  From  this  general  observation  it 
becomes  clear  that  both  the  configuration  of  line-contact  and  two-point-contact 
induce  dryout  over  a  large  area.  However,  limited  dryout  occurs  for  the  configuration 
of  one-point-contact. 

In  practical  heat  transfer  devices  it  is  difficult  to  avoid  or  to  control  the  way  of 
contact  between  the  tube  and  the  support  plate.  But  it  is  possible  to  drill  the  hole 
on  support  plate  with  a  special  configuration ,  for  example,  a  square  hole  that  the 
contact  occurs  between  the  tube  and  the  support-plate  but  with  a  large  opening  for 
flow  by-pass.  To  understand  the  dryout  at  this  condition,  calculations  are  performed 
for  the  geometry  of  line-contact  but  with  various  averaged-gap-thickness  of  the 
annulus,  c.  As  shown  in  Figure  2  the  larger  the  value  of  c  the  more  the  difference  of 
the  radii  of  curvatures  near  the  contact  location.  The  calculated  results  are  shown  in 
Figure  8  which  indicates  that  the  larger  the  difference  of  the  radii  of  curvatures  the 
less  the  dryout  zone.  In  the  present  study,  when  c  is  larger  than  2.117  mm  the 
dryout  zone  will  disappear  completely.  In  fact,  the  similar  conclusion  has  been 
obtained  from  the  experimental  studies  reported  in  the  reference  [1],  A  very  good 
confirmation  between  the  present  analysis  and  the  experimental  results  in  [1] 
suggests  that  this  computation  method  can  be  used  as  an  effective  tool  for  better 
design  of  steam  generator  support-plates  to  minimize  or  avoid  the  dryout,  induced 
corrosion. 
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7.  The  calculated  dryout  pattern  in  annular  clearances  with  inclined  tub* 
contact  at  one-pomt. 
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CONCLUSION 

The  present  subchannel  analysis  gives  reasonable  prediction  of  boiling  and  dryout 
in  confined  space  at  forced  convection.  For  heated  tube  passing  through  round  holes 
on  support-plate  with  small  clearances,  the  dryout  always  occurs  near  the  location  of 
contact.  When  the  difference  in  radii  of  curvature  increases,  the  dryout  zone 
decreases  and  finally  disappears.  Utilizing  this  particular  feature  which  is  predicted 
in  the  present  analysis,  better  design  of  the  support-plate  can  be  achieved. 
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CHAPTER  4 

EXPERIMENTS  OF  BOILING  AND 
DRYOUT  IN  ANNULAR  CREVICES 


WITH  CLOSED  BOTTOM 


INTRODUCTION 


I.  Literature  Review 

Katto  and  Yokoga  et.  al  [1  ]  [2]  [3]  [A]  investigated  confined  space 
boiling  of  saturated  water  on  a  horizontal  flat  plate.  A  movable  optical 
assembly  was  placed  above  a  heated  surface  and  had  a  small  gap  in  between. 
The  boiling  in  the  narrow  gap  has  higher  heat  transfer  coefficient  which  is  due 
to  the  existence  of  thin  liquid  film  under  bubbles  .  They  suggested  that  the 
dryout  condition  in  confined  space  is  a  balance  between  the  consumption  of 
the  liquid  on  the  heated  surfac  and  the  supplying  of  the  liquid  between  the 
intermittent  jetting  of  vapor  mess. 

Ishibashi  and  Nishkawa  [  5]  conducted  cue  boiling  experiment  for  vertical 
annulus  with  open  ends.  It  was  determined  that  there  are  two  distinct  boiling 
regimes,  each  having  different  heat  transfer  characteristics.  No  quantitative 
information  for  dryout  had  been  obtained. 

The  experiment  of  boiling  of  an  horizontal  tube  within  constrain  was 
conducted  by  Jensen,  Cooper  &  Bergers  [  6  ].  The  increase  in  the  heat 
transfer  was  also  explained  by  the  thin  film  evaporation.  The  dryout  heat  flux 
was  found  to  be  directly  proportional  to  the  clearance  and  inversely 
proportional  to  the  length  of  annulus.  Compared  to  a  open  tube,  the  dryout 
heat  flux  was  reduced  by  a  factor  of  10. 


Kusuda  and  imura  [7],  [  g  ]  performed  experiments  with  subcooled  and 
saturated  liquid  for  boiling  heat  transfer  in  vertical  closed  bottom  tubes  and 
annuli.  The  annuli  experiments  consisted  of  an  unheated  core  and  an  externally 
heated  tube.  The  heat  transfer  coefficient  of  boiling  was  reported  to  be 
independent  of  the  length  and  the  diameter  of  the  heated  tube.  The  correlation 
for  dryout  heat  flux  was  attempted  by  dryout  analysis  but  it  was  not  able  to 
correlate  satisfactorily  with  their  experimental  results.  The  visualization  of 
boiling  was  done  in  separate  tests. 

2  Objective 

Increased  knowledge  about  boiling  heat  transfer  in  a  closed  bottom 
annulus  will  improve  our  understanding  of  thermal-hydraulic  related  corrosion 
in  nuclear  steam  generators.  The  available  data  in  the  literature  on  saturated 
boiling  in  vertical  closed  bottom  annulus  is  not  sufficient  for  a  thorough 
analysis. 

It  has  been  demonstrated  by  Ishibashi  and  Nishkawa  [5],  that  in  a 
confined  space,  heat  transfer  characteristics  are  closely  associated  with  the 
boiling  pattern  in  the  crevice.The  experiment  with  visualization  and  heat 
transfer  measurement  at  the  same  time  will  give  more  detailed  understanding 
of  the  boiling  mechanism  in  the  vertical  closed  bottom  confined  space. 

With  the  geometric  configuration  shown  ptWiOUSty  the  crevice  is  likely  to 
have  permear.ant  nucleation  sites  at  its  bottom  end.  Therefore  the  experiment 
will  be  conducted  in  a  closed-bottom  vertical  annulus  with  existing  nucleation 
sites  at  its  bottom.  The  objectives  of  this  experimental  research  arc  as 
follows: 

I.To  investigate  the  effect  of  gap  size  and  system  pressure  on  saturated 
boiling  in  vertical  closed-bottom  annulus  with  existing  nucleation  sites. 
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2. To  study  the  effect  of  fluid  properties  on  the  boiling  heat  transfer 
using  Freon-' 13.  acetone  and  distilled  demineralized  water. 

3. To  observe  the  boiling  behavior  in  a  confined  space  through  movie  and 
to  propose  correlations  to  interpret  the  experimental  results. 

Since  the  dryout  heat  flux  has  its  practical  importance,  attention  will  be 
paid  to  the  correlation  of  the  dryout  heat  flux  to  fluid  properties  and 


geometry. 
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EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

1  Experimental  Apparatus 

Experiments  have  been  conducted  by  boiling  of  Freon-113,  acetone  ana 
distilled  demineralized  water  on  the  outside  of  an  electrically  heated  tube 
which  was  placed  inside  a  an  hollow  quartz  cylinder  with  a  small  gap  in 
between.  The  annulus  is  closed  at  the  bottom  and  open  at  the  top.  The  test 
assembly  was  immersed  in  a  pool  of  saturated  liquid  as  shown  in  Fig.  l.  An 
enlarged  view  of  temperature  measurement  device  is  shown  in  Fig.  2. 

The  heating  tube  used  in  alt  the  experiments  is  stainless  steel  304 
seamless  tubing  25.4mm  O.D.  with  0.71mm  wall  thickness  and  101.6mm  heated 
length.  The  top  of  heating  tube  was  thermally  shrink-fitted  into  a  copper  tuo3 
with  the  same  outside  diameter  and  1.42mm  wall  thickness.  The  bottom  end  of 
the  tube  is  fitted  into  a  copper  block.  Direct  current  passes  though  the 
heating  tube  from  the  upper  copper  tube  to  the  bottom  copper  block. 

As  shown  in  Fig.  3,  the  hollow  quartz  cylinders  are  95.3mm  long  with 
63.5mrr  O.D.  The  inside  diameter  26.04,  27.00  to  30.58mm  to  form  different 
sizes  of  gap  with  the  tube.  The  inside  diameter  is  enlarged  at  the  top  and 
bottom  to  accommodate  the  positioner.  The  inner  and  outer  surfaces  of 
quartz  tube  was  polished  to  50-80  finish  and  is  clear  enough  to  see  through. 

The  three  gap  sizes  tested  are  0.32,  0.80  and  2.58mm,  with  all  of  them 


) 
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Figure  3:  Detail  Drawing  of  Hollow  Quartz  Cylinder 
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76.2mm  long.  There  is  another  test  section  with  0.32mm  gap  but  only  25.4mm 
long.  To  achieve  a  uniform  gap  thickness,  the  top  and  bottom  positioners 
were  installed  to  maintain  the  relative  concentric  position  of  the  heating  tube 
and  the  outer  quartz.  The  bottom  of  the  crevice  was  sealed  by  a  Teflon  ring 
which  does  not  interact  with  Freon-113. 

The  Freon-113  pool  is  contained  by  a  101.6mm  I.D.  and  457.5mm  length 
Pyrex  glass  pipe.  The  pipe  is  closed  at  both  ends  by  flanges  with  Teflon 
seal,  with  beaded  end.  The  allowable  operation  pressure  is  4.4atm.  A  12.7mm 
I.D.  tubing  was  connected  to  the  top  of  the  tank  to  release  the  Freon  vapor 
generated  by  boiling.  The  vapor  was  then  condensed  in  a  water  cooled 
condenser.  The  consumed  Freon-113  liquid  was  refilled  from  the  bottom  of 
tank  by  siphon  from  an  outside  reservoir.  The  saturation  temperature  of  the 
pool  was  maintained  by  an  immersion  heater,  which  was  adjusted  by  a  Variac. 
A  type  k  ungrounded  thermocouple  was  installed  in  the  pool  to  measure  its 
temperature,  and  the  system  pressure  was  monitored  by  a  pressure  gauge.  For 
pressurized  tests  a  compressed  nitrogen  cylinder  was  connected  to  supply  the 
pressure.  For  low  pressure  tests  a  diffusion  vacuum  pump  was  used. 

The  electric  current  passing  through  test  section  was  measured  by  the 
voltage  drop  across  a  calibrated  shunt  using  a  Hewlett-Packard  3465S  digital 
multi-meter.  The  power  was  provided  by  two  direct  current  welders.  A  Miller 
Direct  Current  Arc  Welder  supplies  the  current  in  the  range  of  20-300  amps.  A 
Westinghouse  WSH  Arc  Power  Supply  Direct  Current  supplies  current  between 
200-1000  amps.  The  current  was  conducted  to  the  test  section  through  a 
25.4mm  O.D.  copper  cable  which  was  clamped  to  the  connectcs  of  test 
section.  A  supporting  structure  was  built  to  support  the  weight  of  cables. 

Inside  the  heating  tube,  the  temperatures  can  be  measured  at  different 
locations.  Two  type  k  stainless  sheathed  miniature  thermocouples  with  the 
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0.81mm  diameter  were  pressed  against  the  inner  wall  by  springs  with  a  force 
of  0.48  newtons,  "fie  springs  and  thermocouples  were  mounted  on  a  micarta 
cylindrical  body  which  blocked  the  tube  to  avoid  natural  convection.  The 
thermocouple  was  calibrated  within  the  setup  by  immersing  the  heating  tube  in 
constant  temperature  open  pool.  The  pool  temperatures  were  determined  by 
maintaining  different  liquids  at  their  saturation  temperatures.  The 
thermocouples  were  connected  to  OMEGA  Trendicator.  The  displayed 
temperature  has  an  accuracy  of  ±0.56  °C. 

2  Experimental  Procedure 

The  same  heating  tube  was  used  for  all  the  tests.  The  surface  was 
polished  before  each  test  with  a  #320  sandpaper  to  remove  any  print,  soot, 
flux  or  other  surface  irregularities  so  that  a  consistent  surface  finish  would  be 
obtains  for  all  tests. 

The  liquid  level  was  maintained  nearly  constant  at  100  mm  above  the  top 
of  the  crevice  during  the  test  by  periodically  adding  make-up  liquid  to  the 
pool.  If  the  liquid  is  only  slightly  above  the  heating  tube,  violent  boiling 
could  uncover  the  tube  and  caused  premature  dryout.  The  test  section  is 
preheated  by  a  heat  flux  of  21  kw/m2  and  preboiled  for  at  least  one  hour.  For 
water  tests  with  the  smallest  gap  size,  gas  bubbles  may  be  trapped  inside  the 
annulus  .  The  preheating  for  water  test  required  much  longer  time  then  the 
other  liquids  to  remove  these  entrapped  air  bubbles. 

The  nucleation  site  at  the  bottom  of  crevice  is  formed  naturally  between 
the  Teflon  seal  and  the  copper  block.  In  steam  generator  the  crevices  at  the 
bottom  of  the  tube  sheet  may  contain  permeanent  nucleation  sites.  Therefore 
the  test  is  performed  with  the  nucleation  sites  at  the  bottom  of  crevices. 
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The  steady-state  data  was  taken  3  minutes  after  the  change  of  heat  flux. 
The  thermocouples  are  located  t.2cm  below  the  top  opening  of  auartz  in  order 
to  detect  the  early  dryout.The  movie  was  taken  with  a  Borax  movie  camera  at 
the  speed  of  64  frames  per  second  with  top  lighting.  Kodak  4712  Tungsten 
light  balanced  color  film  was  used. 

The  location  of  thermocouple  could  be  changed  during  the  test  to  check 
the  axial  variation  of  wall  temperature.  The  wall  heat  flux  was  calculated 
from  the  measured  eletric  current  and  the  known  electric  resistance  of  the 
stainless  steel  tube.  The  inner  wall  temperature  was  measured  by 
thermocouple.  With  Known  inside  temperature  the  outer  surface  temperature 
was  calculated  from  the  steady-state  heat  conduction  equation  in  cylindrical 
coordinates,  with  one-dimensional  heat  flow  and  constant  physical  properties. 


I 
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RESULTS  AND  DISCUSSIONS 

1  The  Boiling  Heat  Transfer  of  Freon- 113 

Both  the  confined  space  and  open  pool  experiments  was  conducted  in  the 
Freon-113  pool.  The  results  are  discussed  in  the  following  sections. 

1.1  Boiling  in  the  Open  Pool 

The  tests  were  done  in  the  Freon-113  pool  without  the  quartz  cylinder  at 
its  outside.  The  heat  is  transferred  by  natural  convection  until  boiling  occurs. 

1.  Natural  convection:  At  low  heat  flux,  natural  convection  is  the 

dominant  heat  transfer  mode.  The  natural  convective  heat  transfer  for  a 
vertical  cylinder  with  large  Prantdl  number  is  reported  by  Eckert  [^]  as 

Nu  «  0.508  (Pr/(0.92*Pr)«Gr  Pr)1'4 

X  X 

where  the  x  is  the  distance  measured  from  the  bottom  of  tube.  The 
experimental  data  are  compared  with  this  equation  in  Fig .  The  error  is 
within  10%. 

2.  Pool  boiling: 

At  a  higher  heat  flux,  bubbles  are  generated  on  the  heating  surface.  The 
heat  flux  increases  drastically  with  the  wall  temperature.  In  this  experiment 
the  relationship  between  the  heat  flux  and  wall  temperature  super  heat  is  in 


the  form  of 
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The  data  of  present  F-113  experiments  are  plotted  in  Fig.(ref.(pool 
boiling))  .  Also  included  in  the  same  Figure  is  the  experimental  data  of  Yilmaz 
and  Westwater  [I0]fcr  a  horizontal  copper  tube  (  20  gauge,  6.4mm  O.D.) 
heated  by  steam.  The  curves  of  these  two  experiments  are  parallel  although 
they  are  in  different  orientation  with  respect  to  gravity.  It  is  important  to 
point  out  that  the  open  pool  results  of  present  experiment  does  not  reach  the 
critical  heat  flux. 

The  data  of  another  F-113  experiment  by  Hesse  [11]  also  plotted  in  the 

Fig.  4  .  The  tube  is  made  of  nickel  with  14mm  O.D.  and  heated  by  hot 

water  situated  horizontally  in  the  pool.  The  boiling  curve  is  to  some  extent 
different  from  the  result  of  present  study  and  the  data  of  Yilmaz  and 

Westwater.  The  discrepency  may  due  to  the  different  liquid-soiid  combination 
and  orientation  of  the  heater. 

1.2  Boiling  in  the  Confined  Space 

The  boiling  of  F-113  in  confined  space  has  been  studied  with  the  gap 

sizes  of  0.32,  0.80  and  2.58mm.  The  effect  of  different  lengths  and  operating 
pressures  are  also  studied  for  the  smallest  gap  size.  The  test  conditions  are 
as  listed  in  Table  1. 


Table  1. 

Summary  of  F-113 

Experiments 

Test 

Gap 

Length 

Pressure 

Sections 

(mm) 

(mm) 

(atm) 

1 

0.32 

25.4 

1 

2 

0.32 

76.2 

0.6,1,2.02 

3.04,4.4 

3 

0.80 

76.2 

1 
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4  2.58  76.2  *  1  . 

open 

pool  -  101.6  1 

Boiling  in  the  confined  space  has  been  closely  observed  in  the  present 

study.  The  most  important  feature  is  that  bubbles  are  squeezed  due  to  the 

small  gap  between  the  walls.  The  shape  and  motion  of  the  bubble  are 

seriously  affected  by  the  restricted  space.  In  this  experiment  the  bubbles  are 

mostly  generated  from  the  bottom  of  the  crevice.  With  the  addition  of  heat 

flux  the  phenomenon  of  boiling  in  smalt  gaps,  go  through  the  following  stages: 

F'9-  5. 

1.  Isolated  Deformed  Bubbles: 

At  low  heat  flux  the  bubbles  are  relatively  isolated.  The  shape  of  the 
bubble  is  generally  hemispherical  with  a  flat  bottom.  This  is  in  a  similar 

shape  to  the  very  large  bubble  in  a  open  pool  because  the  drag  effect  is 
important.  A  thin  film  of  liquid  is  formed  between  the  bubble  and  the  walls. 
With  the  addition  of  heat, the  bubble  enlarges  itself  while  moves  up.  At  low 
heat  flux,  a  small  fraction  of  the  heated  area  is  covered  by  bubbles  and  at  a 
higher  heat  flux,  more  area  wili  be  covered. 

2.  Violent  Evaporation  of  Deformed  Bubbles: 

At  a  higher  heat  flux,  the  wall  temperature  super-heat  starts  to  increase. 
The  upward  moving  bubbles  expand  rapidly.  Most  of  the  surface  area  is 

covered  by  the  up-flowmg  deformed  bubbles  with  the  liquid  flowing  down 
from  the  top  in  between  the  bubbles.  The  bubble  rises  with  a  stop  and  go 
pattern  by  merging  or  interacting  with  other  bubbles.  Temporary  dryout  may 
occur  under  the  bubbles  but  liquid  wili  soon  rewet  the  surface.  For  the  low 
temperature  superheat  there  is  no  bubble  nucleation  on  the  solid  surface. 
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3.8ubble  nucleated  on  solid  surface 

If  the  gap  is  large  and  with  a  high  heat  flux,  bubbles  are  generated  in 
the  thin  film  under  the  deformed  bubbles.  The  small  bubble  starts  at  the 
nucleation  site  on  the  solid  surface  and  breaks  through  the  thin  liquid  film. 
Affected  by  the  motion  of  the  large  deformed  bubbles,  the  small  bubble  is 
nucleated  at  random  location.  The  bubble  is  also  generated  outside  the 

area  covered  by  the  large  deformed  bubbles. 

4.  Dryout  Condition: 

Dryout  in  this  experiment  is  defined  as  the  condition  when  the  wall 
temperature  keep  increasing  with  time.  When  the  heat  flux  is  close  to  the 
dryout  value,  the  liquid  flows  downward  as  a  film  at  both  sides  of  the  gap 
the  gas  flows  upward  in  between.  The  liquid  film  on  the  heating  surface  is 
evaporated  extensively  but  the  liquid  film  on  the  unheated  quartz  cylinder  may 
extend  to  the  bottom  of  crevice.  The  flowing  film  at  the  heated  wall  may  dry 
up.  However,  the  liquid  film  along  the  unheated  wall  can  .sometimes,  be 
transferred  t^  the  heated  side  either  by  the  surface  wave  which  bridges  the 
gap  or  pushed  by  the  bubbles  generated  in  the  liquid  pool  which  is 
accumulated  at  the  bottom  of  crevice.  Therefore  the  heated  surface  may  also 
get  cooled. 

The  dryout  heat  flux  for  the  crevice  is  generally  much  lower  than  the 
critical  heat  flux  of  the  open  pool.  The  dryout  heat  flux  also  varies  with  the 
gap  size  and  lengtn.  It  appears  that  the  limiting  criterion  of  dryout  heat  flux 
is  no  longer  the  same  as  that  in  pool  boiling.  The  particular  flow  pattern 
caused  by  the  confinement  effect  has  to  be  considered  as  well. 


5.  Post-Dryout  Condition: 
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Beyond  the  dryout  heat  flux  the  temperature  of  test  section  rises  aiong 
the  time.  To  avoid  the  burn  out  of  the  test  section,  the  power  must  be  shut 
off  in  a  short  duration  of  time.  However,  it  is  possible  tv  atudy  the  post¬ 
dryout  condition  for  crevices  with  very  small  gaps  because  the  dryout  heat 
fluxes  of  these  crevices  are  low  and,  therefore,  the  temperature  excursion  rate 
is  low. 

At  this  condition  the  bulk  of  liquid  at  the  top  of  the  crevice  can  hardly 
get  into  the  crevice.  Standing  waves  of  the  interface  are  observed  near  the 
top  of  crevice.  However,  occasionally  the  wave  breaks  up  and  some  drops 
are  ejected  into  the  crevice.  The  drops  evaporate  while  descending  to  the 
bottom.  At  the  bottom  there  may  exist  a  shallow  liquid  slug  where 
evaporation  and  boiling  may  occur.  All  the  outgoing  vapor  provides  enough 
flow  to  maintain  the  unstable  wave  at  the  top  of  crevice. 

1.3  Boiling  Heat  Transfer 

The  boiling  curves  of  Freon-113  are  shown  in  Fig.  6  .  It  can  be  seen 
that  with  the  increase  of  the  gap  size  the  boiling  curve  of  the  confined  space 
approaches  the  open  pool  boiling  curve.  Boiling  in  the  confined  space  is 
characterized  by  the  deformation  of  bubbles.  The  deformation  of  bubbles  can 
be  characterized  by  the  ratio  of  gap  thickness  and  the  bubble  departure 
diameter,  which  is  generally  known  as  the  Bond  number,  which  may  be  defined 
as 

Bo  =S/<<r/g(/>(-/>B»1/2 

The  Bond  number  encountered  in  the  r-113  experiments  are  summarized 
in  Table  2. 


Table  2  Bond  Numbers  in  F-113  Experiments (1  atm) 


(Kw/m 
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Figure  6.  Boiling  Curve  of  Freon-113  in  the  Closed  Bottom  Annulus 


-  73  - 


Gap (mm) 

Bo 

0.32 

1.56 

0.80 

4.19 

2.58 

13.57 

Prom  the  table,  it  is  reasonable  to  expect  that  the  boiling  curves  may 
differ  much  from  the  open  pool  results  if  the  gap  thickness  is  so  small  that 
the  corresponding  Bond  number  is  in  the  ord^r  of  unity.  Comparing  the  three 
boiling  curves  in  Figure  6  we  may  temporarily  conclude  that  if  the  Bond 
number  is  small  the  effect  of  confinement  would  be  apparent  on  the  boiling 
curve.  Following  this,  the  discussion  of  heat  transfer  may  be  directed  toward 
the  following  three  different  conditions  and  each  of  them  correspond  to  a 
boiling  phenomenoa  as  descibed  in  the  last  section: 

(1)  Large  Bond  Number  and  high  heat  flux.(Bubbles  Nucleated  on  Solid 
Surface) 

(2)  Small  Bond  Number  and  low  heat  flux.dsolated  Deformed  Bubbles) 

(3)  Small  Bond  Number  and  high  heat  flux.lViolent  Evaporation  of 
Deformed  Bubbles) 

The  particular  features  of  heat  transfer  at  each  condition  will  be 

discussed  in  turn. 

(1)  Large  Bond  Number 

When  the  Bond  Number  is  large  the  boiling  curve  is  similar  to  that  ot  the 
open  pool.  In  general,  the  wall  super  heat  is  reduced  comparing  to  pool 

boiling.  At  high  heat  flux,  the  boiling  curve  is  almost  in  parallel  to  the  open 

pool  boiling  curve  which  is  in  the  form  of 
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Chernobyl'skii  and  Tananaiko  [U]  performed  experiments  using  water  in 
an  annulus  with  open  ends  .  The  dimension  of  the  annulus  was  varied  to  give 
a  range  of  Bond  number  from  7.9  to  31.7.  They  observed  a  similar  behavior 
of  boiling  curve  as  indicated  in  Figure  6  for  the  large  crevice  which  has  a 
Bond  number  13.57.  They  suggested  that  the  increase  of  heat  transfer 
coefficient  in  the  narrow  space  is  due  to  the  reduction  of  the  bubble  size  and 
the  increase  of  turbulence  in  the  mixture.  The  liquid  is  violently  agitated  so 
that  a  higher  rate  of  heat  removal  from  the  heated  wall  occurs. 

The  same  explanation  can  be  applied  to  the  case  of  the  vertical  annulus 
with  closed  bottom  which  was  studied  in  this  experiment. 

(2)  Small  Bond  Number  and  Low  heat  flux  condition: 

At  this  condition,  the  boiling  phenomena  is  mainly  isolated  deformed 
bubbles.  The  heat  flux  increases  with  little  increase  of  wall  temperature.  The 
wall  temperature  super  heat  is  in  the  order  of  0.56±°C. 

4 

The  temperature  super  heat  of  bubble  can  be  calculated  follow  the 
procedure  of  Griffith  and  Wallis  ref.  [ll].  A  small  pressure  difference  is 
required  to  enable  the  bubDle  to  exist  in  liquid.  The  Clausius-Clapeyron 
equation  is  used  to  calculate  the  bubble  super  heat  from  the  pressure 
difference.  In  this  experiment  bubble  is  flattened  with  a  radius  of  curvature 
approximately  half  of  the  gap  size.  The  bubble  super  heat  becomes  as 

AT  *  2T  <t/H  p  S 
sat  fgrv 

The  calculated  superheated  temperature  of  a  flat  F-113  bubble  is  0.005°  C 
Therefore  the  low  superheated  temperature  in  the  deformed 


bubble  is  possible. 
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In  this  low  heat  flux  the  deformed  bubble  only  covers  certain  percentage 
of  heated  area. If  alt  the  heat  is  used  to  evaporate  the  liquid  into  the  bubble 
,the  total  heat  flux  should  equal  to  the  heat  flux  under  the  deformed  bubble 
but  corrected  by  the  time-averaged  area  ratio  (Ar)  of  deformed  bubbles  on 
heated  surface. 


q  *  Arq 
w  Mb 


The  heat  flux  under  a  deformed  bubble  can  be  evaluated  by  the  increases 
of  bubble  size  over  a  period  of  time.  Assuming  all  the  heat  transfer  under 
the  bubble  is  used  to  evaporate  the  liquid,  the  heat  flux  under  the  deformed 
bubble  can  be  shown  as 


q  *  AV/>  H  /AM 
Mb  '  v  tg 


The  heat  flux  under  the  bubble,  the  velocity  of  bubble  and  the  area  ratio 
are  summarized  in  Table  3. 


Table  3.  Sur.n-.ary  of  Results  from  Visual  observation 


Total  Heat 

Time  Interval  between 

Time  Averaged 

Heat  Flux  j 

Flux 

Successive  Figures 

Area  Ratio 

under  Bubble 

q  „ 

At 

Ar 

q  i 

D 

(kw^m2) 

(sec.) 

(kw/m2) 

0.62 

0.09 

32% 

2.88 

0.62 

0.03 

32% 

0.81  !  ! 

0.62 

0.125 

32% 

1.33  i 

C.62 

0.125 

32% 

i 

0.21 

0.16 

16% 

1.29 

0.21 

0.125 

16% 

1.41 

I 
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0.21 


0.125 


16% 


The  total  heat  flux  equals  to  the  thin  film  evaporative  heat  flux  times  the 
area  ratio  covered  by  the  thin  film.  With  the  area  ratio  and  the  average  heat 
flux  measured  from  figures  the  total  heat  flux  calculated  are  0.54  kw/m2  and 
0.22  kw/m2  which  compared  with  0.62  kw/m2  and  0.21  kw/m2  obtained  from 
electric  power  the  error  is  within  ±13%. 

From  the  experimental  data  it  can  be  concluded  that  most  of  the  heat  is 
used  in  the  evaporation  of  the  deformed  bubbles.  The  variation  of  the  total 
heat  flux  changes  the  area  covered  by  deformed  bubbles, but  not  changing  the 
wall  superheat  temperature.  The  boiling  curve  becomes  almost  a  vertical  line 
on  the  heat  flux  versus  wall  temperature  .  That  is  what  we  observed  on  the 
boiling  curve,  (Fig.  £  )  under  the  condition  of  small  Bond  number  and  low  heat 
flux. 


3)  Small  Bond  Number  and  High  heat  flux: 

At  this  condition  most  of  the  heated  surface  is  covered  by  deformed 
bubbles.  Unstable  dryout  under  bubble  is  observed  in  high  heat  flux.  The  wall 
temperature  increases  with  the  increase  of  heat  flux.  The  relationship  is  direct 
proportional 

q*AT 

At  a  fixed  wall  temperature  the  heat  flux  is  higher  with  smaller  gap  size. 
If  the  Bond  number  is  smaller  than  1.84,  this  region  may  not  be  observed 
easily. 

At  high  heat  flux  the  coalition  of  buoble  occurs  frequently.  Eventually  a 
center  vapor  core  is  formed  with  wavy  liquid  film  flowing  down  on  both  side. 


77  - 


2  Boiling  Heat  transfer  of  Acetone  and  Water 

In  additional  to  Freon-113,  boiling  in  acetone  and  distilled  demmeralizd 
water  is  studied.  The  major  effect  in  boiling  heat  transfer  using  different 
fluids  is  the  change  of  physical  properties.  Since  the  Bond  number  appears  as 
an  important  parameter  in  describing  the  heat  transfer  in  confined  space,  the 
corresponding  Bond  numbers  of  each  fluid  in  various  crevices  are  listed  in  the 
Table  4. 

The  acetone  and  water  experiments  extend  the  range  of  the  Bond  number 
to  a  lower  range  than  that  involved  in  Freon-113  experiments.  Therefore  the 
boiling  curves  of  acetone  and  water  may  be  more  in  similar  to  a  that  of  small 
crevice  in  Freon-113  experiments. 


Table  4.  Summary  of  Bond  number  of  present  test 


Gap/llcucd 

(nun) 

F-113 

Acetone 

water 

0.32 

1.66 

1.03 

0.72 

0.80 

4.19 

2.60 

1.84 

2.56 

13.57 

8.23 

5.87 

The  boiling  curves  of  acetone  are  shown  in  the  Figure  y  .  For  crevices 
with  Bond  number  less  than  8.23  the  curves  are  similar  to  those  of  small  Bond 
numbers  in  Freon-113  experiments.  The  boiling  curve  of  the  large  gap  has  a 
Bond  number  8.23  and  behaves  similar  to  F-113  large  gap  curve  with  Bo  equal 
to  13.57  and  quite  different  from  the  others. 

For  the  water  experiments  all  the  conditions  have  the  Bond  number  less 
or  equal  to  5.87.  Therefore  the  effect  of  confinement  would  be  significant. 
As  shown  in  the  Figure 


I 
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all  the  boiling  curves  in  the  present  study  have  similar  characteristic  to 
the  curves  of  low  Bond  numcer  m  Freon-113  or  acetone  experiments. 


The  relationship  among  the  Bond  number,  the  heat  transfer  characteristics, 
and  the  bubble  region  are  shown  in  the  following: 


Table  5.  The  Bond  Mo.  with  respect  to  Bubble  Region  and 
Associated  Heat  Transfer  Characteristics 

Bond  Number 

Bubble  region  0.72-1.66  1.84-5.78  8.4  13.6  Infinite 


Isolated  Deformed 

Bubbles (Small  5cnd  X  XX 

Me.  Low  Heat  Flux) 

Violent  Evaporation 

of  Deformed  Bubbles  X  XX 

(Small  Bond  No.  High 
Heat  Flux) 

Bubble  Generated  under 

Deformed  Bubbles (Large  X  X 

Bond  Number) 

Open  Pool  X 


From  Table  5  the  Bond  number  is  a  parameter  to  indicate  the  impor.-.nce 
of  boiling  regions  for  a  given  system.  The  larger  the  Bond  number  ;s.  the 
closer  it  is  to  open  pool  results. 

3  The  Dryout  Heat  Flux 

As  shown  in  previous  sections  the  dryout  heat  flux  in  the  crevice  with 
the  closed  bottom  is  considerably  lower  than  the  critical  heat  flux  for  boiling 
in  a  open  pool.  For  Freon-113  experiment,  the  critical  heat  flux  of  a  open 
pool  boiling  is  201  kw/m2,  but  for  the  smallest  gap  (0.32mm)  the  dryout  heat 
flux  studied  presently  is  only  6.3  kw/m2  ,  which  is  almost  2  order  of 
magnitude  less.  The  value  of  dryout  heat  flux  is  summarized  in  Table  6. 
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Generally,  the  smaller  the  gap  size  the  lower  the  dryout  heat  flux.  The  wall 
super  heat  at  dryout  13  also  fo«nc  to  be  decreasing  with  smaller  gap  size. 
The  dryout  heat  flux  decreases  when  the  length  of  the  crevice  increases. 

Table  6.  The  Summary  of  Dryout  Heat  Flux  Data  in  Crevices 


S(mm) 

L(ram) 

P (atm) 

q  (Kw/m2) 

0 

0 

Water 

2.58 

76.2 

1 

210.1 

16.11 

0.80 

It 

N 

82.0 

9.72 

0.32 

M 

It 

25.2 

3.06 

Acetone 

2.58 

76.2 

It 

109.5 

7.78 

0.8C 

It 

It 

37.2 

7.78 

0.32 

25.4 

It 

35.0 

18.33 

M 

76.2 

It 

10.8 

3.05 

F- 112 

2.58 

76.2 

It 

6S.4 

18.06 

0.80 

•t 

m 

21.1 

15.83 

0.32 

25.4 

•t 

17.4 

4.72 

n 

76.2 

w 

6.29 

0.83 

M 

* 

0.60 

5.54 

1.11 

M 

tt 

2.02 

7.55 

0.63 

m 

tt 

3.04 

8.53 

1.39 

m 

•» 

4.04 

CO 

t- 

• 

00 

1.67 

1.  Analysis 

The  significant  difference  between  the  pool  boiling  critical  heat  flux  and 
the  dryout  heat  flux  of  crevices  may  attributed  to  the  different  limiting 
mechanism  for  the  liquid  to  come  into  the  heating  surface.  Near  the  dryout, 
the  outward  vapor  flow  prohibit  sufficient  liquid  to  flow  down  into  the  crevice 
Fig.  9  _  This  phenomena  is  similar  to  that  of  hydrodynamic  counter-current 
flooding  in  tubes.  When  flooding  occurs  the  down-coming  liquid  becomes 
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stagnant  because  the  gravity  force  is  balanced  by  the  up  going  vapor  shear 
force. 


The  counter-current  flooding  of  two  phase  flow  in  a  round  tube  is  a  well 
studied  problem.  For  tubes  with  small  diameter  the  correlation  proposed  by 
Wallis  [  14]  base  on  the  superficial  velocities  is  the  best  fit  of  data,  the 
hydrodynamic  flooding  experiment  of  the  annulus  geometry  with  Bond  number 
equals  25  has  been  conducted  by  Shire  et  al  [Iff].  The  results  can  also  be 
correlated  satisfactorily  by  a  similar  form  of  Wallis  correlation  but  with 
different  empirical  constant. 

On  the  other  hand,  the  dryout  heat  flux  for  tubes  of  different  diameters 
has  been  correlated  by  Nejat  using  flooding  considerations.  In  the  derivation, 
Nejat  [16]  used  the  Wallis  correlation  together  with  the  limiting  condition  for 
large  tubes.  He  argued  that  the  dryout  heat  flux  for  tubes  with  small  diameter 
is  a  function  of  Bond  number,  and  the  correlated  experimental  data  have  the 
Bond  number  in  the  range  of  18  to  100. 

In  the  presen*  experiment  the  range  of  Bond  number  for  Freon-113, 
acetone,  and  distilled,  demineralized  water  is  between  0.72  and  13.6.  Since  the 
Bond  number  is  generally  small,  the  Wallis  correlation  may  be  used  to 
describe  the  hydrodynamic  counter-current  flooding  in  crevice.  For  an  annulus 
with  closed  bottom  the  flooding  will  start  at  the  top  opening  where  the  vapor 
velocity  is  the  highest.  At  a  higher  heat  flux  the  dryout  starts  at  the  top  of 
crevice.  It  is  also  important  to  notice  that  the  dryout  condition  may  occur 
before  the  flooding  happens  because  the  dryout  occurs  when  the  supply  of 
liquid  is  "insufficient"  to  remove  the  heat  applied  to  the  tube.  At  this 
moment  the  liquid  flow  has  not  been  completely  stopped  yet. 


The  Wallis  correlation  can  be  represented  as 
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j  *u2+  j/,/2-  C  if  Bo  <  50 

where;  )a*-  Vfl,/2^D,/>«Vfl^1/2 

V*  W'^v^r1'2 

and 

i  *  M  lip  A) 

g  9  9 

j,  *  M}/(/j(A) 

where  the  characteristic  dimension  D,  which  is  measured  from  the  center 
of  the  gap,  is  used  to  non-dimensionalize  the  superficial  velocities.  Since  the 
annulus  has  a  closed  bottom,  the  steady  state  mass  flow  rate  of  liquid  come 
in  at  the  top  opening  of  the  crevice  should  equal  to  the  mass  flow  rate  of 
vapor  flowing  out  of  the  top  opening. 

M 

9  ' 

At  the  moment  the  dryout  starts,  the  rate  liquid  flows  into  the  crev.ce 
equal  the  rate  liquid  evaporates.  The  mass  velocity  can  be  related  to  the  heat 
flux  as 

Mg  ■  q  A*/H 
w  fg 

By  substituting  and  rearranging  the  equation^  and  assuming 

the  gap  size  is  much  smaller  than  the  diameter,  the  dryout  heat  flux  can  be 
represented  as 

B,x(M/>  Ip  ),/4)2x.  Ip*  1C2 
l  ri  rg  rg  rt 

where  B  *q  L/H  p  VS 

I  \v  tsr9 
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and  V  -p  ~1/2(gD(/oiv,»1/2 

'  3 

The  boiling  number  Bo  is  defined  as  a  ratio  of  outgoing  vapor  velocity 
(q  L/H  />  S)  to  the  bubble  rising  velocity(/>  _,/2(gD(/>  -/>  ))1/2)  in  the  annulus 

w  fgr  g  t  '  8 

[  1 4]  .  The  characteristic  dimension  for  the  annulus  is  the  diameter  of  the 
annulus  D  which  is  measured  at  the  middle  of  the  gap.  In  terms  of  the 
diameter  of  the  inner  tube  D  and  the  gap  thickness  S  the  D  is 

D  *  Di  *  S 

The  dryout  heat  flux  of  present  study  for  Freon-113,  acetone  and  distilled 
demineralized  water  in  crevices  of  different  gaps  and  heights  can  be  correlated 
by  the  above  formulation.  In  addition,  the  water  experimental  data  of  Kusuda 
&  Imura  [8  3  with  Bond  number  3.17-18  is  also  correlated.  The  results  are 
shown  in  Figure  10. 

All  these  data  can  be  correlated  properly  in  the  form  of 

B  x(1»(/>  ip  )W4)2xo  Ip  -  0  38 
i  r\  q  g  • 

It  is  important  to  include  the  results  of  Kusuda  &  Imura  into  the 
correlation.  In  their  experiment  the  annulus  is  heated  from  outside  wall. 
However,  the  present  experiments  a^e  heated  from  inside  tube.  The  dryout  heat 
flux  correlation  is  independent  of  the  side  which  is  heated. 

The  crevice  of  the  smallest  size  (0.32mm)  has  also  been  tested  in  Freon- 
113  pool  at  different  pressures.  In  general,  the  behavior  of  boiling  at  other 
pressures  is  similar  to  that  discussed  previously.  At  elevated  pressure  the 
bubbles  appear  to  be  slightly  more  toward  spherical  shape.  The  dryout  heat 
flux  is  also  fitted  well  by  the  proposed  correlation. 


It  is  generally  observed  that  the  Wallis  hydrodynamic  flooding  correlation 
for  small  gaps  is  independent  of  the  surface  tension.  Therefore,  the  present 
dryout  heat  flux  correlation  is,  similarly,  independent  of  the  surface  tension. 
Although  the  water  has  3  times  larger  surface  tension  than  that  of  Freon-113 
and  acetone  these  data  can  be  correlated  satisfactorily  by  the  equation  which 
does  not  containing  the  surface  tension. 


CONCLUSIONS  AND  RECOMMENDATIONS 


1  Conclusions 

1.  From  the  visual  observation  the  boiling  in  the  confined  space  has  three 
different  boiling  regions:  isolated  deformed  bubble,  violent  evaporation  of 
deformed  bubble  and  bubble  nucleated  on  the  solid  surface.  Bond  number  can 
be  used  as  a  parameter  to  indicate  the  existence  of  boiling  regions  for  a 
giving  system. 

2.  For  the  small  Bond  number  and  low  heat  flux,  the  heat  flux  increases 
with  little  increase  of  wall  temperature.  The  boiling  is  characterized  as  many 
isolated  deformed  bubbles.  At  higher  heat  flux,  the  wall  temperature  increases 
proportional  ’o  heat  flux  and  the  boiling  is  under  violent  evaporation  of 
deformed  bubble.  For  large  Bond  number,  heat  flux  can  go  even  higher  with 
small  bubbles  with  bubble  nucleated  on  the  solid  surface.  The  heat  transfer 
characteristic  is  similar  to  pool  boiling  curve  but  with  slight  less  wall 
superheat.  The  Larger  the  Bond  number  is.  The  closer  the  results  to  pool 
boiling  curve. 

3.  The  heat  transfer  under  the  isolated  deformed  bubble  is  dominated  by 

thin  film  evaporation.  The  increasing  of  heat  flux  increases  the  area  covered  by 
deformed  bubbles, but  not  changing  the  wall  superheat  temperature. 

4.  Under  the  condition  of  having  bubbles  nucleated  on  solid 
surface  _  the  recjuctjon  0f  .wa||  superheat  comparing  to  open  pool  boiling  curve 
is  attributed  to  the  convective  effect. 
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5.  The  dryout  heat  flux  in.  the  confined  space  is  much  lower  than  that  in 
the  pool  boiling.  The  correlation  for  dryout  heat  flux,  derived  from  the 
counter-current  flooding  condition,  fit  the  experimental  data  within  ±24  %  for 
various  liquids  at  different  conditions.  For  annulus  heated  from  inside  tube  or 
outside  wall  the  dryout  heat  flux  can  be  correlated  with  the  same  correlation 
proposed  here. 

2  Recommendations  for  Further  Study 

1.  Boiling  with  isolated  deformed  bubble  shows  very  good  heat  transfer 
characteristics.  The  further  study  on  this  phenomena  may  result  in  better 
understanding  of  heat  transfer  augmentation  in  confined  space. 

2.  For  the  purpose  of  understanding  the  thermohydraulic  induced  corrosion 
in  the  heat  exchanger  of  PWR,  the  dryout  heat  flux  for  eccentric  annulus 
should  be  studied.  It  is  the  likely  case  in  the  real  heat  exchangers. 

3.  The  dryout  heat  flux  of  confined  space  with  both  ends  is  suspected  to 
be  a  phenomenon  of  counter-current  flooding  also.  For  open  ends  the  liquid 
may  enter  the  confined  space  form  the  bottom  end.  The  conservation  of 
mass  used  in  our  analysis  of  closed  bottom  annulus  may  not  be  applied  directly. 
With  some  modification  in  the  derivation  of  equations  and  carefully  perform 
experiments  similar  type  of  dryout  correlation  may  be  obtained. 

4.  The  heat  exchanger  in  the  nuclear  submarine  is  placed  horizontally. 
The  same  corrosion  problem  may  occur.  It  is  valuable  to  study  the  horizontal 
closed  bottom  annulus.  It  is  suspected  that  the  boiling  behavior  of  isolated 
deformed  bubble  remains  the  same. 

6.  This  test  was  done  with  an  active  nucleation  site  at  the  bottom,  it 
will  be  interesting  to  study  the  boiling  in  the  confined  space  without  such  a 


sjte  exists. 
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CHAPTER  5 


EXPERIMENTS  OF  CONVECTIVE  BOILING  AND 
DRYOUT  IN  ANNULAR  CREVICES 


Boiling  of  dryout  in  annular  crevices  with  forced  convective  is 
under  study  currently.  The  loop  has  been  constructed.  The  schematic 
of  the  loop  is  shown  in  Fig.  1  with  the  loop  characteristics  summarized 
in  it.  In  order  to  perform  boiling  experiments  at  constant  pressures, 
an  accumulator  is  attached  to  the  loop  to  absorb  the  volume  expansion  of 
the  two  phase  fluid  duirng  the  tests.  Large  flow  by-pass  is  provided  to 
the  test  section  to  maintain  the  flow  stability  in  the  test  section. 

The  instrumentation  system  has  also  been  established. 

The  schematic  of  the  test  section  is  shown  in  Fig.  2.  The  quartz 
test  section  is  changeable  for  different  lengths  and  annular  gaps. 

The  data  of  the  smallest  gap  annulus  are  shown  in  Fig.  3.  The  general 
feature  is  similar  to  that  of  conventional  boiling  curve  of  large-gap 
annulus  but  with  relatively  low  wall  superheat  and  low  critical  heat  flux. 
The  trend  of  the  critical  heat  flux  with  respect  to  the  mass  flux 

also  similar  to  the  conventional  boiling  however 
with  lower  values  of  critical  heat  flux. 

Other  data  are  undertaking  and  more  detailed  results  and  discussion 


will  be  available  soon. 
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